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NASA  Kennedy  Space  Center  (KSC)  carried  out  an  analysis  of  the  effects  of  aeroacoustics 
produced  by  stationary  solid  rocket  motors  in  processing  areas  at  KSC.  In  the  current 
paper,  attention  is  directed  toward  the  acoustic  effects  of  a  motor  burning  within  the  Vehicle 
Assembly  Building  (VAB).  The  analysis  was  carried  out  with  support  from  ASRC  Aerospace 
who  modeled  transmission  effects  into  surrounding  facilities.  Calculations  were  done  using 
semi-analytical  models  for  both  aeroacoustics  and  transmission.  From  the  results  it  was 
concluded  that  acoustic  hazards  in  proximity  to  the  source  of  ignition  and  plume  can  be 
severe;  acoustic  hazards  in  the  far-field  are  significantly  lower. 


Nomenclature 

A  area,  m2 

Aw  A-weighted  sound  pressure  level  correction,  dB 

a  speed  of  sound,  m/s 

b  line  support  spacing,  m 

C  model  coefficient 

c  phase  velocity,  m/s 

d  diameter,  m  or  interior  wall  spacing,  m 

DI  directivity  index,  dB 

F  thrust  or  force,  N 

/  frequency,  Hz 

L  length,  m 

Lw  overall  sound  pressure  level,  dB 

M  Mach  number 

m  mass  per  unit  area,  kg/m2 

ra  mass  flow  rate,  kg/s 

OASPL  overall  sound  pressure  level,  dB 

P,  p  pressure  or  overpressure,  Pa 

R  gas  constant,  J/kg-K  or  distance  to  receiver  from  plume  source,  m 

Ra  A-weighted  frequency  sensitivity  weighting  function,  W 
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total  absorption  surface  area,  m2 

Strouhal  Number 

sound  pressure  level,  dB 

plume  path  length  from  nozzle  to  source,  m 

temperature,  K 

transmission  loss,  dB 

thickness,  m 

velocity,  m/s 

sound  power,  W  or  weight,  lb 
width,  m 

overall  acoustic  sound  power,  W 
plume  core  length,  m 
absorption  coefficient 

geometric  angle  from  receiver  to  local  plume  slice  flow  direction,  radians 
effective  angle  corrected  for  Strouhal  Number  effects,  radians 
bandwidth  of  frequency  band,  Hz 
plume  slice  length,  m 

acoustic  efficiency,  decimal  or  wall/ window-panel  material  loss  factor 

specific  heat  ratio 

wavelength,  m 

Poisson's  Ratio 

density,  kg/m3 

modulus  of  elasticity,  psi 

incidence  angle,  radians 

transmissivity 

Subscripts 

ambient 
point  b 
nozzle/bore 
chamber  or  critical 
point  d 

nozzle  exit  or  point  e 
frequency  or  point  f 
incident 

outer  sheet  of  sandwich 
plume  slice 

plume  slice,  frequency 

reflected 

total 

nozzle/bore,  blast-origin,  or  fundamental 
lower  critical 
higher  critical 
sonic  condition 


I.  Introduction 

r  J^HE  Vehicle  Assembly  Building  (VAB)  was  built  to  allow  processing  of  Apollo  Program  launch  vehicles  in  the 
1960’s  and  later  was  adapted  to  process  the  Space  Shuttle  launch  vehicles.  Since  then  the  VAB  has  served 
effectively  as  the  location  to  assemble  and  mate  solid  rocket  boosters  (SRB)  to  the  vehicles  External  Tank  (ET)  and 
Orbiter.  The  VAB  will  again  be  adapted  to  process  new  human  and  non-human  rated  launch  vehicles  in  support  of 


NASA’s  strategic  plan  for  space  exploration. 


Previously,  for  the  Space  Shuttle  Program  (SSP),  the  Kennedy  Space  Center  (KSC)  Safety  Office  determined  a  safe 
separation  distance  (SSD)  of  1,315  feet  for  the  inhabited  building  distance  (IBD)  [7].  The  calculation  used  the 
Department  of  Defense  (DoD)  explosive  safety  standards  (sometimes  referred  to  as  the  weight  based  approach).  The 
SSP  safe  separation  distance  for  the  VAB,  derived  using  the  DoD  weight  based  approach,  allows  for  processing  of 
up  to  four  complete  Shuttle  SRB’s  containing  approximately  4.44  million  pounds  of  solid  propellant.  The  currently 
proposed  Constellation  Program  (CxP)  launch  schedule  could  require  assembly  and  storage  of  nearly  13  million 
pounds  of  hazard  classification  1.3  propellants;  up  to  eight  5-segment  ARES-V  boosters  to  be  stored  in  the  VAB 
simultaneously.  Using  the  weight  based  approach  for  eight  ARES-V  boosters  would  result  in  a  safe  separation 
distance  of  1,810  ft.  This  distance  extends  beyond  a  number  of  high  occupancy  facilities  in  the  area  around  the 
VAB.  The  NESC  concluded  that  the  DoD  weight  based  approach  was  overly  conservative  for  SRMs  with  a  1.3 
hazard  classification.  NASA  Engineering  and  Safety  Center  (NESC)  recommended  that  an  alternative  approach 
based  on  calculating  the  threat  due  to  each  hazard  component  (heat  flux,  toxics,  acoustics,  etc.)  be  pursued  [1]. 

In  the  event  of  an  inadvertent  ignition,  several  hazardous  components  could  either  prevent  or  restrict  egress  of 
personnel  and  potentially  impact  the  surrounding  facilities.  A  methodology  development  as  well  as  an  assessment  of 
the  noise  levels  in  proximity  of  the  VAB  was  required  to  support  the  alternative  approach  recommended  by  the 
NESC.  Acoustic  SSD’s  as  well  as  levels  of  resistance  nearby  facilities  have  from  acoustic  energy  was  addressed. 
Several  scenarios  were  considered  and  later  in  the  study  the  case  set  was  bounded  using  results  from  parallel  studies. 
These  parallel  studies  helped  narrow  the  VAB  building  configuration  and  ignition  propagation  as  a  result  of  an 
inadvertent  ignition.  The  emphasis  of  this  study  was  placed  on  the  far-field  sound  hazards  to  personnel  and  the 
effectiveness  of  surrounding  facilities  as  safe  heavens  from  acoustic  energy. 

Prediction  of  aero-acoustic  noise  as  well  as  transmission  into  nearby  facilities  was  performed  using  semi-analytical 
models.  The  facilities  modeled  were  those  in  immediate  proximity  to  the  VAB;  the  Launch  Control  Center  (LCC) 
and  Orbiter  Processing  Facility  #3  (OPF-3).  An  exposure  criterion  for  personnel  was  developed  to  establish  safe 
separation  distances.  All  methodologies  and  findings  were  examined  by  a  NASA  peer  review  committee. 


II.  Modeling  Assumptions 


Currently,  CxP  would  reuse  the  VAB  with  minor  modifications  to  support  processing  of  the  Ares-I/V  vehicles.  The 
VAB  is  the  fourth  largest  building  in  the  world  by  volume.  It  consists  of  4  High  Bays  used  to  process  NASA  launch 
vehicles,  each  assembled  on  a  Mobile  Launch  Platform  (MLP).  Several  platform  levels  are  used  to  perform 
processing  operations  and  will  too  be  reused  and/or  modified.  Each  HB  is  separated  by  concrete  walls  (lightly 
reinforced)  and  large  steel  columns  which  provide  the  main  support  to  the  structure. 


Figure  1:  NASA  Ares-I/V  Flight  Vehicles  and  the  VAB 

The  shell  is  made  up  of  hundreds  of  aluminum  “punch-out”  panels  which  are  designed  to  fail  above  approximately 
100  psf  (meant  to  protect  the  facility  from  large  pressure  loading  during  severe  weather).  It  was  shown  from  CFD 
models  that  the  pressure  loading  from  ignition  over-pressure  (IOP)  and  plume  impingement  (resulting  from  an 
inadvertent  ignition)  would  knock  over  internal  walls  and  fail  all  VAB  paneling.  In  summary,  shortly  after  ignition 
only  the  VAB  ceiling,  main  supports,  and  ML’s  would  remain  intact;  it  was  assumed  the  vehicles  remain  intact.  This 
assertion  was  used  in  the  acoustic  modeling  and  helped  to  simplify  the  analysis  such  to  make  the  acoustic  fields 
diffuse  into  a  fourth  space. 

Though  a  fully  processed  vehicle  would  certainly  become  propulsive  if  an  inadvertent  ignition  occurred,  an  un¬ 
capped  booster  would  not.  This  is  because  the  chamber  pressure  isn’t  sufficient  enough  to  create  enough  thrust.  To 
simplify  the  analysis,  it  was  assumed  that  all  burning  boosters  were  fixed  in  place  in  the  VAB.  Furthermore  nozzle 
gimbling  and  motor  pitching  was  ignored. 


Furthermore  it  was  assumed  that  plumes  emanating  from  boosters  and  impinging  on  the  VAB  floor/ceiling  spread 
symmetrically.  This  assumption  made  the  acoustic  calculations  easier  because  the  fields  could  be  treated  as  two- 
dimensional.  Also,  in  scenarios  were  adjacent  boosters  (on  the  same  MLP)  burned  together  (and  in  the  same 
direction),  their  plume  acoustic  powers  were  combined  and  an  equivalent  plume  was  computed.  This  simplification 
is  typically  done  in  aeroacoustics  since  the  plumes  become  conjoined  if  in  close  proximity. 


Sound  absorption  (dissipation)  from  the  atmosphere  is  dependent  on  temperature  and  relative  humidity.  The 


atmosphere  at  KSC  is  fairly  humid  year  round  (typically  40-95%)  and  temperature  can  vary  between  extremes 


(typically  35-105°  F).  Above  approximately  4  kHz,  temperature  and  relative  humidity  become  competing  effects. 


Atmospheric  Absorption  (Orlando  Florida) 


Figure  2:  Atmospheric  Sound  Absorption  versus  Relative  Humidity  (for  various  frequencies  and  temperatures) 


The  dissipation  from  atmospheric  absorption  can  be  significant  above  5  kHz.  For  this  study  the  absorption  was 


ignored  to  be  conservative  and  simplify  the  analysis.  Furthermore,  ground  reflections  other  than  the  initial  plume 


deflection  were  ignored. 


III.  Theoretical  Background 


The  computational  modeling  of  aeroacoustics  from  plumes  is  difficult  because  of  the  small  time  steps  and  large  flow 
times  required  to  resolve  spectra  (10  kHz  ~  10"4  seconds).  Furthermore,  models  used  to  convert  the  unsteady- 


turbulent  fields  to  acoustic  sources  are  computationally  expensive.  Alternatively,  semi-empirical  methods  have  been 
used  that  combine  theoretical  equations  with  empirical  correlations.  The  method  used  in  this  study  breaks  the  plume 
into  several  slices  which  are  converted  to  monopole  sources  (see  figure  below). 


Source  Locations 


Figure  3:  Plume  conversion  into  acoustic  sources  [2] 

The  sources  are  then  characterized  and  sound  is  radiated  inside  the  domain.  This  method  has 
agreement  with  measured  data.  Most  recent,  comparisons  with  data  from  the  Ares-IX  launch 
agreement  (see  a  summary  below). 


Proximity 

Spectrum  Error 

OASPL  Error 

Near-field;  ML/Launch  Pad  (0-500’) 

Mid-field;  Launch  Site  (1000-1400’) 

±  0.5  dB 

+  1  dB 

Far-field;  adjacent  pad/VAB  (5000-10,000’) 

±5  dB 

Figure  4:  Model  error  determined  from  Ares-IX  launch  test  data 


IV.  Exposure  Criteria 

Several  sources  stating  the  sound  pressure  levels  associated  with  human  exposure  exist.  OSHA  standard  1019.95  [3] 
provides  exposure  limits  for  industrial  areas  based  on  exposure  time.  Because  capped  and  uncapped  boosters  can 
bum  for  long  periods  of  time  (>  5  minutes),  using  the  OSHA  table  would  result  in  a  115  dBA  exposure  limit.  The 
intent  of  the  OSHA  table  was  for  use  in  industrial  environments  and  not  for  inadvertent  ignition-type  events. 
Furthermore,  use  of  the  OSHA  value  would  result  in  a  large  SSD  (overly  conservative).  Hearing  damage  will  largely 
affect  a  person’s  ability  to  egress  because  of  the  ear's  sensitivity. 


proven  to  give  good 
showed  phenomenal 


A-Weighted  Correction 


Figure  5:  A-Weighted  Correction  [4] 

Ruptured  ear  drum(s)  will  impede  a  person’s  balance  and  ability  to  egress;  this  occurs  around  160  dB  [4].  The  ear’s 
threshold  of  pain  occurs  around  140  dB;  intense  nausea  is  typical  at  levels  around  and  above  158  dB  [4]. 

The  personnel  exposure  limit  used  in  this  study  was  140  dB  (corresponding  to  the  threshold  of  pain).  At  these  levels 
personnel  will  experience  a  large  degree  of  permanent  and  irreversible  hearing  damage  but  not  ear  drum  rupture; 
thus  allowing  the  ability  to  egress.  Because  of  the  log(R 2)  nature  of  acoustic  sound  propagation,  revising  the 
criteria  upward  will  drastically  reduce  the  SSD.  Also,  OASPL  levels  are  measured  at  a  head  level  of  5 ’6”  and 
revision  of  this  target  height  to  a  larger  value  would  also  drastically  reduce  the  SSD. 

V.  Calculation  Methodology 

Acoustic  predictions  were  performed  using  a  modified  Eldred’s  2nd  Method.  The  plume  is  broken  into  several  slices 
and  given  a  unique  spectrum.  The  sound  pressure  from  all  the  slices  (at  each  octave  band  frequency)  is  summed  to 
determine  the  net  sound  pressure  at  a  specific  point.  This  is  done  first  by  solving  for  the  plume  metrics  using  the 
chamber  conditions.  Using  provided  chamber  conditions  and  nozzle  geometry,  nozzle  exit  conditions  were  found 
assuming  ID  compressible  isentropic  flow  [6].  Assuming  the  flow  to  be  choked  at  the  nozzle  throat,  mass  flow  rate 


is  found, 


m  = 


pcA*  In  /  2  x(rc+i)/(rc-i) 

^  \yc  +  l) 


The  exit  Mach  number  is  found  using  the  area-Mach  number  relation, 


(yc+i)/(7c-i) 


Eq.  1 


Eq.  2 


The  equation  has  two  solutions;  subsonic  and  supersonic  condition.  Since  the  chamber  pressure  is  such  that  it’s 
above  the  critical  pressure  ratio,  we  assume  the  supersonic  solution. 


Area-Mach  Number  Relation 


Figure  6:  Area-Mach  number  relation  (for  various  specific  heat  ratios) 
The  exit  plane  pressure  is  then  found  using  the  Mach  number, 


-yc/(yc-i) 


Knowing  the  exit  Mach  number,  the  exit  plane  temperature  is  then  found, 


=r«( 


With  the  exit  temperature  found,  exit  plane  speed  of  sound  is  found, 

& e  sfcKJe 


Eq.  3 


Eq.  4 


Eq.  5 


The  exit  plane  velocity  is  found  then  using  the  Mach  number  and  speed  of  sound, 


Thrust  is  then  found  to  be  [9], 


F  =  mVe  +  ( Pe  -  Pamb)Ae  Eq.  7 

Plume  core  length  is  found  using  the  exit  Mach  number-diameter  correlation  [2], 

=  3.45(1  +  0.38 Me)2  Eq.  8 

dg 

This  correlation  was  derived  from  various  experimental  data  sets  of  various  plume  types;  it  was  a  least-square -type 
fit  to  that  data. 


Figure  7:  Non-dimensional  plume  core  length 

In  order  to  calculate  the  overall  sound  power,  an  acoustic  efficiency  must  be  used.  With  direct  floor  impingement, 
the  efficiency  is  found  using  the  data  in  SP-8072  in  conjunction  with  the  nozzle  geometry  and  placement  relative  to 
the  VAB  floor  and/or  ceiling. 
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Figure  8:  Acoustic  efficiency  of  various  plume  deflections  [2] 

The  overall  sound  power  is  then  found  using  the  efficiency  factored  against  the  mechanical  power  [4], 


W0A  =  VFVe  Eq.  9 

The  overall  sound  power  can  then  be  translated  into  an  overall  sound  pressure  level  [2], 

Lw  =  10  log(WOA)  +  120  Eq.  10 

Once  the  plume  metrics  are  known,  the  plume  is  discretized  into  several  slices  and  Strouhal  numbers  are  found  for 
each  plume  slice  using  the  flow  variables  and  choosing  a  frequency  range  (typically  10°-104  Hz), 

f  S  CL  e 


Sn(s)  = 

ve  U-amb 

The  normalize  overall  sound  power  of  each  slice  is  found, 


Eq.  11 


loiog^^Vi^J  =  wl°g 


C4\C  5 


Eq.  12 


where  Ci,  c2,  c3,  c4,  and  c5  are  correlation  coefficients  (omitted  for  SBU  considerations).  The  strongest  source 
typically  sits  between  1.5-2  core  lengths  from  the  exit  plane  (shown  in  figure  below).  The  core  (inviscid)  region 
terminates  (xt)  when  the  jet-edge  shear  layers  collapse  and  form  a  large  turbulent  zone  in  the  plume.  This  is  also  the 
cause  of  the  log  reduction  in  axial  velocity  along  the  plume  axis. 


Normalized  Acoustic  Power  in  Plume 


Figure  9:  Normalized  overall  sound  power  along  plume 
The  overall  acoustic  power  for  each  slice  is  found  [2], 

xtW(sy 


L w,ps  10  log 
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/A  xt(s)\ 
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Knowing  the  Strouhal  number  range,  the  normalized  power  spectrum  is  found  for  each  plume  slice, 

10109  -  moS  (T  cfcfs%,) 


Eq.  13 


Eq.  14 


where  c1?  c2,  c3,  and  c4  are  correlation  coefficients  (omitted  for  SBU  considerations).  An  important  caveat  to  these 
empirical  correlations  is  that  their  coefficients  are  dependent  on  the  propulsion  type  (chemical,  nuclear,  etc). 


Plume  Slice  Normalized  Acoustic  Power  Spectrum 


Figure  10:  Normalized  power  spectrum  of plume  slice 
The  power  spectrum  for  each  plume  slice  is  then  found  [2], 

lw,ps,f  =  10  log  (w(f,s)/W{s)f^})  +  Lwps  -  10  log  +  i0lOg  (A/)  Eq.  15 

\  sae  /  \ s aamb / 

The  directivity  indices  are  found  by  determining  the  true  angles  the  receiver  location  makes  with  each  plume  slice. 
The  true  angles  are  then  converted  to  equivalent  angles  correcting  for  Strouhal  number  effects, 


Peff  =  P  ~  Cilog  (c2Sn(s))  Eq.  16 

where  Ci  and  c2  are  correlation  coefficients  (omitted  for  SBU  considerations).  Knowing  the  effective  angle, 
directivity  index  is  computed, 


D1{Peff)  =  10  log 


_ Cl  +  C2COS4(jSe//) _ 

[(1  -  c3cos(/?e//))c«  +  C5]e6[l  +  c7exp(cspeff)\ 


-  c9log{Sn(s))  -  cw 


Eq.  17 


where  c1?  c2,  c3,  c4,  c5,  c6,  c7,  c8,  c9,  and  Ci0  are  correlation  coefficients  (omitted  for  SBU  considerations).  The 
average  (over  a  typical  Strouhal  number  range)  directivity  of  chemical  rockets  it’s  approximately  40°  (see  figure 
below). 


Plume  Source  Directivity  Index 
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Figure  11:  Acoustic  source  directivity  index  (at  various  Strouhal  numbers) 

Again,  the  coefficients  vary  for  different  propulsion  types.  The  mean  directivity  for  various  propulsion  sources  is 
shown  in  the  figure  below. 


Figure  12:  Directional  characteristics  for  various  propulsion  types  [2] 


The  frequency  dependant  sound  pressure  level  from  each  plume  slice  for  a  receiver  is  [2], 


SPLps,f  =  Lwpsf  -  10  log(i?2)  -  11  +  Dl(peff) 


Eq.  18 


A  logarithmic  summation  over  all  plume  slices  of  the  above  will  yield  the  sound  pressure  spectrum  for  a  receiver 

[2], 

/  SPLpSj\ 

SPLf  =  10  log  I  )  10  io  j  Eq.  19 


The  overall  sound  pressure  level  for  a  receiver  can  be  found  from  another  logarithmic  summation  over  the  entire 
spectrum  [2], 


OASPL  = 


Eq.  20 


Corrections  to  the  sound  pressure  spectrum  for  ear  sensitivity  can  be  done  by  including  the  A-weighted  term  with 
the  SPLps  f  equation, 


SPLf,Aw  =  SPLf  +  Aw(f) 


The  A-weighted  correction  for  the  human  ear  is  frequency  dependant  [4], 


Ra(D 


122002/4 

(/2  +  20.62)V(/2  +  107 .72)(/2  +  737.9 2)(/2  +  122002) 


KiD  =  2-0  +  20  log(RA(n) 


Eq.  21 


Eq.  22 
Eq.  23 


Summation  of  SPLf, aw  will  give  the  OASPL  in  dBA. 


Corrections  to  the  sound  pressure  spectrum  (normal  not  A-weighted)  for  diffraction  with  a  wall  can  be  done  by 
including  the  diffraction  term  with  the  SPLpS;f  equation.  An  empirical  correlation  used  to  determine  the  diffracted 
gain  in  the  acoustic  spectrum  (see  figure  below). 
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Figure  13:  Diffracted  gain  versus  non-dimensional  characteristic  length  [5] 


The  non-dimensional  characteristic  length  can  be  converted  to  frequency  by, 


Eq.  24 


The  acoustic  spectrums  are  then  corrected  for  diffraction  by  applying  the  correlation, 


SPLftdiff=SPLf  +  DIFF{n 


Eq.  25 


The  diffracted  acoustic  spectrums  are  then  applied  to  all  walls  of  the  facility. 


SPLf, inner  =  SPLf,diff  ~  10  loglO  (=)  +  10  logi 


i4(l  -  a ) 


Eq.  26 


In  the  above  equation,  the  second  term  is  the  transmission  loss  (TL;  from  wall  and  windows)  and  third  is  the 
absorption  loss  (from  internal  carpeting,  cubicles,  etc). 


The  average  transmissivity  is  found  by  taking  an  area-weighted  average  for  all  walls.  The  transmission  loss  is 
depended  on  the  wall  fundamental  and  critical  frequency, 
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Eq.  27 
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Eq.  28 


The  transmission  loss  is  also  largely  dependent  on  the  wall  type  (single  or  double  panel  wall).  For  a  single  panel 


wall, 
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It  is  important  to  note  that  the  method  conservatively  assumes  there  is  no  transmission  loss  below  the  fundamental 
frequency  of  the  panel.  This  was  done  primarily  because  specific  treatment  of  this  region  was  not  taken  into  account 
being  that  the  sound  field  is  diffuse  [11]. 
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Figure  14:  Typical  transmission  loss  profile  for  a  single  panel  wall 
For  double  wall  constructions,  the  analysis  methodology  is  significantly  different.  First,  the  calculations  for  the 
method  must  be  done  in  metric  units.  This  is  due  to  several  simplifications  that  were  included  in  the  development  of 


the  procedure  by  Bies  and  Hansen  [11].  Secondly,  since  double  wall  constructions  can  contain  an  air  space,  the 
equation  for  the  fundamental  frequency  changes  to  account  for  coupling  effects  [8], 


/o  =  80 
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Eq.  30 


d(m1m2) 

Once  this  frequency  is  determined,  it  then  becomes  necessary  to  determine  the  critical  frequency  of  both  wall  panels 

[8], 


L  2  (12  (l  -  (yjfjj  gq  3] 

where  if  any  wall  panel  is  made  of  thin,  equal-thickness  sheets  of  material  sandwiching  a  lightweight  core,  the 
critical  frequency  equation  changes  to  [8], 

c,  L(2(l-(.,4  Eq.  32 

2?rj  El,2tl,2(tPl2-tl,2) 

It  is  important  to  note  that  the  subscripts  1  and  2  in  these  equations  correspond  to  the  panel  with  the  lowest  critical 
frequency  and  highest  critical  frequency,  respectively.  The  general  transmission  loss  spectrum  for  a  double  panel 
construction  is  shown  in  the  figure  below. 
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Figure  15:  Typical  transmission  loss  profile  for  double  panel  wall 


For  a  double  panel  wall  the  transmission  loss  is  found  from, 
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Eq.  33 


-78  dB 


0.5 fc2  </<fc2 


LINEAR  INTERPOLATION  USED 


Cl=B  +  6  + 1  -  log10  (tj2  )  dB  (if  fc2  *  fcl  ) 

C2=B  +  6  + 1  -  log10  (rj2 )+  51og10  (rjx )  dB  (if  fc2  =  fcl ) 


fc2  <  f  —ft 


LINEAR  INTERPOLATION  USED 


F  =  - 


15 


login  (2) 


(log10  (//  )-  log10  (fc2 ))+  c  dB 


ff  >fc  2 


An  important  aspect  of  the  double  wall  calculation  method  is  that  linear  interpolation  is  required  to  solve  for  the 
transmission  loss  at  the  frequency  values  between  the  known  points.  This  is  due  to  the  fact  that  there  are  several 
possible  combinations  of  equations  that  need  to  be  used,  depending  on  the  boundary  conditions  of  the  wall/window 
being  analyzed.  An  example  of  this  is  when  calculating  the  transmission  loss  at  point  B  on  the  spectrum.  If  the  space 
between  the  two  panels  comprising  the  double  wall  construction  contains  no  sound  absorbing  material,  the  value  of 
B  i  is  used  at  point  B;  otherwise,  the  value  at  point  B  is  the  larger  of  B !  and  B2. 


The  average  absorption  coefficient  is  found  by  taking  an  area- weighted  average  of  all  individual  absorption  sources, 


a  = 


aiAi+ai+1Ai+1+--+anAn 

Ai+Ai+1+--+An 


Eq.  34 


Typical  values  for  absorption  coefficient  are  shown  in  the  table  below. 


ACOUSTIC  ABSORPTION  COEFFICIENTS 

MATERIALS 

125  HZ 

250  HZ 

500  HZ 

1000  HZ 

2000  HZ 

4000  HZ 

Painted  Concrete: 

0.1 

0.05 

0.06 

0.07 

0.09 

0.08 

Carpet: 

0.01 

0.05 

0.10 

0.20 

0.45 

0.65 

Window  Glass: 

0.18 

0.06 

0.04 

0.03 

0.02 

0.02 

Linoleum  Floor: 

0.02 

0.03 

0.03 

0.03 

0.03 

0.02 

Ceiling  Tiles: 

0.30 

0.51 

0.80 

0.85 

0.78 

0.66 

Gypsum  Board 

0.29 

0.10 

0.06 

0.05 

0.04 

0.04 

Removable  Walls 

0.14 

0.10 

0.06 

0.04 

0.04 

0.03 

Figure  16:  Typical  values  of  absorption  coefficient  [10] 


VI.  Sample  Scenarios 

Since  CxP  is  the  current  NASA  program,  a  sample  case  for  the  larger  Ares-V  heavy  lift  vehicle  follows.  In  the 
scenario,  an  Ares-V  vehicle  is  in  the  final  stages  (capping)  of  SRB  stacking  when  an  inadvertent  ignition  occurs. 
This  event  could  happen  in  any  of  the  four  HB’s.  Because  of  this,  the  results  are  superimposed  on  all  HB  locations 
to  determine  the  final  SSD. 


Figure  17:  Sample  case  booster/plume  configuration 


The  SRB  segments  are  assumed  to  be  identical  to  the  SSP  segments  for  this  sample  case.  It’s  assumed  that  the 
chamber  condition  of  these  incomplete  stacks  is  approximately  100  psia  and  5000°  F;  the  chamber  gas  specific  heat 
ratio  was  assumed  1.17.  Nozzle  geometry  used  in  the  SSP  SRB’s  was  assumed  along  with  the  same  SRB/ET  sitting 
position  on  the  ML  in  the  VAB. 


Because  the  boosters  are  uncapped,  two  plumes  develop.  The  lower  plume,  because  of  the  nozzle,  has  a  larger  exit 
Mach  number  (Me  ~  2)  plume  than  the  top  plume  (Me  =  1).  As  a  result,  the  lower  plume  sources  are  distributed 
farther  out  into  the  farfield  and  radiate  sound  further.  In  the  near-field,  the  upper  plume  sound  radiation  is  mainly 
confined  to  the  VAB;  little  of  its  sound  propagates  to  any  significant  area  outside  the  VAB  (see  the  figure  below). 


OASPL  Contour  (dB)  for  VAB  Nearfield 
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Figure  18:  Near-field  OASPL  contour  for  sample  case 

In  the  farfield,  the  sound  radiation  is  dominated  mainly  by  the  lower  plume.  After  approximately  500  feet  the 
OASPL  decreases  rapidly  until  it  drops  to  insignificant  levels  past  about  1000  feet  (see  the  figure  below). 
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OASPL  Contour  (dB)  for  VAB  Farfield  (truncated  to  0-1 00’  elevation) 
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Figure  19:  Farfield  OASPL  contour  for  sample  case 

The  SSD  is  drawn  below  with  the  SSP  QD  arc  discussed  earlier.  The  SSD  in  the  sample  case  sits  well  within  the 
previous  SSP  QD  arc  (see  the  figure  below). 


Figure  20:  VAB  area  SSD  plot  for  sample  case  (OASPLLCc  =  148  dB,  OASPL0pf-3  =  145  dB) 

From  the  transmission  loss  analysis,  the  LCC  and  OPF-3  impedances  were  not  enough  to  reduce  the  OASPL  below 
the  140  dB  personnel  limit;  as  a  result  these  facilities  would  be  unsafe  as  havens. 


VII.  Comments  and  Conclusions 

A  methodology  for  determining  a  SSD  for  aeroacoustic  hazards  has  been  successfully  developed.  Future 
improvements  of  the  method  could  aim  to  reduce  the  conservatism  (likely  in  the  range  of  10-20  dB).  A  new  module 
is  planned  such  to  model  the  effectiveness  of  ear  protection  equipment.  Furthermore,  a  module  to  parametrically 
model  the  effectiveness  of  facility  barrier  shielding  is  likely.  The  addition  of  these  two  new  modules  would  add  the 
dimension  of  mitigation  to  the  model  and  methodology. 
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Introduction 

•  In  early  2007  efforts  initiated  to  quantify  inadvertent  hazards  in  KSC 
processing  areas  for  Constellation  Program  (CxP) 

•  Use  of  the  weight-based  approach  methodology  resulted  in  very  large  arcs 

•  Alternatively,  hazards  were  laid  out  individually  to  be  assessed 

•  Some  of  the  hazards  considered  were  blast/fragmentation,  radiative 
heating,  aeroacoustics,  toxics 

•  End  result  was  to  determine  a  safe  separation  distance  (SSD)  for  each 
individual  hazard 

•  The  inputs,  assumptions,  criterion,  and  results  of  each  individual  hazard 
analysis  were  reviewed  by  peer  committee 

•  This  presentation  deals  specifically  with  the  Aeroacoustic  hazard 
component 


ENGINEERING _ 

Assumptions 

•  Burning  stack/motor(s)  are  fixed  in  place  and  not  pitched/tilted 

•  Assumed  VAB  paneling  and  internal-dividing  walls  removed 

-  simplifies  analysis 

-  provides  most  conservative  acoustic  fields  for  farfield 

•  Plume(s)  treated  as  axisymmetric 

-  asymmetric  spreading  of  the  plume  is  not  considered 

•  Adjacent  plumes  in  close  proximity  conjoined 


-  only  if  plumes  within  1-5  exit  diameters 

-  sound  powers  combined  to  create  one  plume 
•  Neglected  atmospheric  sound  absorption 

-  effect  negligible  for  small  regions  considered 

-  provides  small  degree  of  conservatism  to  analysis 


Neglected  ground  reflections  (other  than  initial  plume  deflection) 

Assumed  incident  acoustic  load  on  facility  normal  and  symmetric 
Facilities  simplified  to  rectangles  for  diffraction  analysis 
Walls  and  windows  treated  as  simply  supported  panels 
At  frequencies  where  absorption  data  unavailable,  assumed  no  absorption 


-  difficult  to  find  building  material  data  below  125  Hz 

-  majority  of  rocket  noise  sits  in  the  101  ~  103  Hz  area 

-  provides  a  large  level  of  conservatism  in  analysis 
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Model  Overview 


Source  Locations 
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Each  source  has  spectrum  dependant  on  plume  strength  (in  form  on  Strouhal  number) 

Normalized  converted  to  overall  using, 

tw.psJ  =  10 log(w(f,  s)/W(s )  [^22*] )  +  LWips  -  10 log  (Ys^SS±)  +  iOlog  (A f) 

Directivity  of  source  with  receiver  dependant  on  Strouhal  number  and  true  angle  with  plume  flow 
direction;  SPL  at  frequency  for  given  receiver  and  source, 

M w  =  ~  10log(«J)  -  11  + 


For  any  given  receiver,  sound  pressure  from  each  source  summed  at  each  frequency;  resulting  in  SPL 
spectrum  for  receiver, 

,  /V 

SPLf  =  IQlog  (  y  10  1 


Receiver  SPL  spectrum  integrated  to  find  OASPL 

/V- '  SPLf 

OASPL  -  IQlog  (  y  10  io  J 
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The  diffraction-gain  spectrum  is  computed  from  semi-empirical  correlations  using  the  building 

x  =  108.354  m 


0.] 


Lc/X  - 


1 .0 
M  I,  n  Smurri 
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OPF  3 


OPF  3 


->  x  =87.180  m 

x  =208.957  m 

x  =254.139  m 


WoveJtngih 

The  spectrums  at  the  facility  distances  are  then  corrected  using  the  diffraction-gain  spectrum  to  find 
the  acoustic  excitation  on  the  facility  wall  for  the  transmission  analysis 


SPL\f)  =  SPL(f)  +  A 

diffraction  (/) 

Sample  Acoustic  Spectrum 
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Based  on  the  facility  geometry  and  composition,  wall 
transmission  and  internal  absorption  losses  are  calculated 

/1\  ..  ( A(l  —  a) 

TL  =  10  log10  (=j 


AL  =  10  log 


10 


Sa 


Internal  spectrum  then  computed  by  taking  external  diffracted 
spectrums  and  correcting  for  transmission/absorption  losses 

SPLf  =  SPL*  -TL-\-AL 


r 

Diffracted  Spectrum 


ACOUSTIC  ADSORPTION  COEFFICIENTS 

MATERIALS 

125  HZ 

250  HZ 

500  HZ 

1000  HZ 

2000  HZ 

4000  HZ 

F'aifitctl  £  Irnirrefe; 

tU 

0.05 

0.06 

a.07 

0.05 

D.UH 

Carpet: 

0.01 

0.05 

010 

0.20 

0.45 

0.65 

Window  Glass: 

cus 

0.06 

004 

0.03 

0.02 

0.03 

Linoleum  Hooi: 

0.02 

0.03 

0.03 

0.03 

0.03 

0.02 

Ceiling  Tiles: 

0.30 

0.51 

0.S0 

0.65 

0.78 

0.66 

Gypsum  litr.Trd 

0  20 

0.10 

006 

0.05 

0.04 

0.04 

Rjcitiovablc  Walls 

0  L4 

0.10 

006 

0.04 

0.04 

0.03 

TL  (dB) 


V 


Y 


/.  /,  0.5/,.  f,2  f, 
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Aeroacoustic  Validation 


•  Modified  Eldred's  2nd  Method  used  for  analysis 

•  Original  method  (NASA  SP-8072)  has  been  updated  with  new  core 
and  directivity  correlations 

•  Super-position  principle  has  been  used  to  incorporate  multiple  non- 
conjoined  plumes 

•  Method  used  to  predict  acoustics  around  launch  pad  under  STS 
launches 

•  Method  also  used  to  predict  launch  acoustics  for  Ares-I/X 

•  Sample  data  shown  to  right  for  Ares-IX  test  flight  on  October  28, 
2009  (data  values  not  shown  for  ITAR/SBU  considerations) 

•  Excellent  agreement  between  model  and  measurements  in  the  near 
and  midfield  spectrums;  ±  0.5  dB 

•  Farfield  spectrums  from  model  over  predict  slightly  in  range  of 
interest  (10-1000  Hz);  ±  7  dB 

•  OASPL  values  (not  shown  for  ITAR/SBU  considerations)  between 
model  and  measurements  for  nearfield,  midfield,  and  farfield  off  by 
only  +  1  dB 
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Personnel  Exposure  Limit 

•  OSHA  1910.95(b)(2)  would  require  a  115  dBA  limit 

—  OSHA  table  meant  for  industrial  working  environments  and  not  inadvertent  ignition  events 
-  Would  result  in  overly  conservative  safe  separation  distances 

•  Practical  limit  would  be  that  associated  with  onset  of  physical  impediment 

•  Typical  human  ear  threshold  of  pain  occurs  at  140  dB;  ear  drum  rupture  is  typical  at  160 
dB 

•  Personnel  exposure  limit  used  was  140  dB  (OASPL);  taken  at  conservative  head  level  (5') 


OSHA  1910.95(b)(2) 


Duration  Per  Day,  In  Hours 

Sound  Level  dBA,  Slow  Response 

a 

90 

6 

92 

4 

95 

3 

97 

2 

100 

1-1/2 

102 

1 

105 

1/2 

110 

1/4  or  less 

115 

Note:  When  the  daily  noise  exposure  is  composed  of  two  or  more  periods  of  noise  exposure  of 
different  levels,  their  combined  effect  should  be  considered,  rather  than  the  individual  effect  of  each. 
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Calculation  Methodology 
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Diffraction/Transmission  Procedure 


Sample  Scenario 


•  VAB  is  4th  largest  building  in  the  world  by  volume 

•  Consists  of  4  High  Bays  (HB)  used  to  process  vehicles 
(Shuttle,  Ares-I,  Ares-V) 

•  Ares-V  vehicle  uses  (2)  SRBs  strapped  to  center  stage 

•  SRBs  are  assembled  first  in  VAB  on  ML  piece-by-piece 

•  Sample  scenario  considers  inadvertent  ignition  of  (2) 
4.5-segment  stacks  in  a  single  HB 

•  Event  could  occur  in  any  one  of  the  HBs,  so  SSD 
calculated  and  placed  at  each  HB  center 

•  Final  arc  drawn  from  centroid  of  HBs  such  that  all 
SSDs  encompassed 
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OPF-3  OASPL  =  145  dB  (UNSAFE) 


RED  ARC  =  sample  case  GREEN  ARC  =  SSP 


OASPL  Contour  (dB)  for  VAB  Nearfield 


500- 


400 


Upper  Plufne  Flow  Direction  (ceiling  deflection) 


:  2ooi 


1001 


Lower  Plufrie  Flow  Direction  (floor  deflection) 


1170 


1 160 


-150 


140 


130 
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LCC  OASPL  =  148  dB  (UNSAFE) 

OASPL  Contour  (dB)  for  VAB  Farfield  (truncated  to  0-100’  elevation) 


0  500  1000  1500 


x-coordinate  (ft) 
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Future  Work 

•  Obtain  low  frequency  (<  125  Hz)  data  for  facility  materials 

-  either  by  experimental  testing  or  computationally  (FEA) 

-  remove  some  conservatism  from  analysis 

-  could  be  as  much  as  2-3  dB  for  certain  facilities 

•  Develop  ear  plug/muff  module  using  ANSI  S12. 68-2007  to  better  determine  mitigation 
levels 

-  better  determine  effectiveness  of  safety  protection 

-  better  quantify  effectiveness  of  certain  equipment  for  rocket  acoustic  spectrums  (10-1000  Hz) 

•  Parametric  analysis  of  facility  shielding 

-  overlapping  mitigation  strategy  from  all  hazard  analyses  involved  blast  shielding  facilities  in  close 
proximity  to  VAB 

-  analysis  would  determine  effectiveness  of  shield  size/construction  to  acoustic  hazards 

•  Develop  GUI  for  aeroacoustics  code  and  register  in  NASA  database 

-  easier  to  obtain  funding  to  support  work  if  registered 

-  could  make  available  to  other  agencies  or  possibly  open-source  (ITAR  restricted) 
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BACKUP  SLIDES 
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General  Acoustics 
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SOUND  POWER  LEVEL 


(W\ 

L-  = 101og  fe) 

where  W  is  the  acoustic  power  (typically  in  aero-acoustics  some  percentage  of  the  mechanical  power), 

and  W0  is  the  reference  power  (Watts) 


SOUND  PRESSURE  LEVEL 

SPL  =  20  log  (— ) 

W 

where  p  is  the  root  mean  square  (rms)  pressure, 
and  p0  is  the  reference  rms  pressure  (Pa) 


p0  is  by  standard  2x1 0-5  Pa,  this  is  equivalent  to  the  approximate  minimum 


Prms  <Pa> 

i both  Lw  and  SPL  are  indicated  in  dB 
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Source  of  sound 

Sound  pressure 

(Pascal) 

Sound  pressure  level 

(dB  re:  20  pPa) 

Theoretical  limit  for  undistorted  sound  at  1  atmosphere 
environmental  pressure 

101,325  Pa 

194.0937  dB 

Krakatoa  explosion  at  100  miles  (160  km)  in  air 

20,000  Pa 

180  dB 

Simple  open-ended  thermoacoustic  device 

12,000  Pa 

176  dB 

Ml  Garand  being  fired  at  1  m 

5,000  Pa 

168  dB 

Jet  engine  at  30  m 

630  Pa 

150  dB 

Rifle  being  fired  at  1  m 

200  Pa 

140  dB 

Threshold  of  pain 

100  Pa 

130  dB 

Hearing  damage  (due  to  short-term  exposure) 

20  Pa 

approx.  120  dB 

Jet  at  100  m 

6  -  200  Pa 

110-  140  dB 

Jack  hammer  at  1  m 

2  Pa 

approx.  100  dB 

Hearing  damage  (due  to  long-term  exposure) 

6*10-1  Pa 

approx.  85  dB 

Major  road  at  1 0  m 

2xio-i -6X10-1  Pa 

80  -  90  dB 

Passenger  car  at  1 0  m 

2X10-2-2X10-1  Pa 

60  -  80  dB 

TV  (set  at  home  level)  at  1  m 

2x10-2  Pa 

approx.  60  dB 

Normal  talking  at  1  m 

2x10-3-2x10-2  Pa 

40  -  60  dB 

Very  calm  room 

2xio-4-6x10-4  Pa 

20  -  30  dB 

Leaves  rustling,  calm  breathing 

6x10-5  Pa 

10  dB 

Auditory  threshold  at  2  kHz 

2x10-5  Pa 

0  dB 
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A-Weighted  Correction 


-  low  frequency  and  high  frequency  sounds  perceived  to  be  not  as  loud  as  mid-frequency  sounds 

-  ear  most  sensitive  to  noise  around  2-6  kHz  range 

-A-weighted  curve  is  the  standard  for  quantifying  sounds  pressure  levels 
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Aero-Acoustics 
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Overall  Sound  Pressure  Level  Map 


Decreasing  Plume  Deflection 
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Normalized  Acoustic  Power  in  Plume 


25 


KSC 

ENGINEERING 


increasing  nozzle  expansion  ratio 
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Plume  Source  Directivity  Index 


1-10 

s 


-15 


0  20  40  60  80  100  120  140  160  180 

Exhuast  Relative  Angle  (deg.) 


Dl{patt)  =  lOlog 


0.37  -f  0.37 


1(1  -  -I-  0.1637] : 'll  + 


-  0.698EoS(STIl»)  -  1.67 


— S  =0.001 

n 

— S  =0.01 

n 

— S  =0.1 

n 

S  =1.0 
n 
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Acoustic  Diffraction 
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Diffraction  Overview 

•  Bending  or  spreading  of  waves  around  obstacle 

-  larger  wavelengths  will  bend  around  smaller  objects  or  barriers,  but  not  vice  versa 

•  Smaller  wavelengths  will  not  diffract  around  smaller  barriers  as  effectively 

-  result  known  as  'shadow  zone'  where  the  sound  is  quieter  than  elsewhere  around  object 

•  Same  conclusion  can  be  made  for  long  wavelengths/longer  barriers,  which  result  in 
a  'shadow  zone'  due  to  reflection 

•  Result  of  diffraction  is  reflection,  which  can  affect  incident  wave  (enhance  sound 
intensity  levels  on  surface) 


Total  Pressure  =  Incident  Pressure  +  Pressure  Reflected 

PTOT  =  Pi  +  Pr 
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Diffraction  Analysis 

Wyle  dB  Loss  Model 


NOTE:  Wyle  model  is  average  diffraction  correction  to  be  applied;  if  the  value  of  Lc  /  A  is 

greater  than  10,  the  dB  gain  is  assumed  to  be  6  dB. 


30 


Ksc 

ENGINEERING 


Acoustic  Transmission 
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EFFECTS  OF  ACOUSTIC  LOADS  ON  STRUCTURES 


Sound  can  either  be  reflected,  absorbed,  or  transmitted  through  a  wall 


ACOUSTIC  ABSORPTION  COEFFICIENTS 

MATERIALS 

125  HZ 

250  HZ 

500  HZ 

1000  HZ 

2000  HZ 

4000  HZ 

5”  Thick  Fiberglass 

0.50 

0.97 

0.99 

0.99 

0.97 

0.93 

Carpet  on  Concrete 

0.11 

0.14 

0.24 

0.40 

0.50 

0.60 

Heavy  Window  Glass 

0.18 

0.06 

0.04 

0.03 

0.02 

0.02 

Marble  Tile 

0.01 

0.01 

0.01 

0.01 

0.02 

0.02 
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SOUND  TRANSMISSION  LOSS 
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•  Transmission  loss  represents  acoustic  resistance  of  structure 

•  Allows  you  to  quantify  the  drop  in  sound  through  a  wall 


splx  • 


TL 


%  SPL2 


SPL2  =SPLl-TL  +  101og10 


Sa  J 


TL  =  -101og10(r) 


SINGLE  WALL  SOUND  TRANSMISSION  LOSS 


DOUBLE  WALL  SOUND  TRANSMISSION  LOSS 


TL  hlB) 


log  ,„(///,) 


log  ,0  (/) 
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